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Abstract 

We consider dijet production at a e^e'^ collider in a class of effective theories with the 
relevant operators being four-fermion contact interaction. Despite the nonrenormalizable 
nature of the interaction, we explicitly demonstrate that calculating QCD corrections is 
both possible and meaningful. Calculating the corrections for various differential distribu- 
tions, we show that these can be substantial and significantly different from those within 
the SM. Furthermore, the corrections have a very distinctive flavor dependence. 
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1 Introduction 

Measurement of multi-jets rates in electron-positron annihilation provide an excellent test of 
perturbative quantum chromodynamics (QCD) [1]. In Standard Model (SM), QCD has been 
tested in the perturbative regime to a high degree of accuracy [2]. Though SM is very success- 
ful model in high energy physics, there are theoretical issues that cannot even be addressed 
within the framework of the SM alone. Examples include the replication of the fermion fami- 
lies, the naturalness problem associated with the Higgs scale, charge quantization, the baryon 
asymmetry in the universe, the presence of dark matter etc.. Clearly, an answer to such vital 
questions may be obtained only in a model much more ambitious than the SM. Candidates for 
the role include, amongst others, supersymmetry [3], grand unification [4,5] (with or without 
supersymmetry) , family symmetries (gauged or otherwise) and compositeness for quarks and 
leptons [12]. Therefore the possible existence of new physics beyond SM, involving the four 
fermion contact interaction would be one of the viable model and may give rise to small effects 
at present and future linear collider. We illustrate this explicitly in the case of dijet production 
via electron-positron annihilation process. We limit our calculation for exclusive two-jet cross 
section in this model. 

The replication of fermion families in SM suggests the possibility of quark-lepton composite 
structure or bound states of more fundamental constituents, called preons [6]. 't Hooft has 
shown that the gauge theory [7] of preon binding naturally reproduce the composite fermions 
of less massive than the preon binding scale which is called characteristic scale A. At this scale 
A, this interaction would become very strong leading to bound states (composites) which are 
to be identified as quarks and leptons. In most such models [8,9], quarks and leptons share at 
least some common constituents. Since the confining force mediates interactions between such 
constituents. Well below the scale A, such interactions would hkely be manifested through an 
effective four fermion contact interaction [10] term that is an invariant under the SM gauge 
group. A convenient and general parametrization of such interactions is given by [11, 12] 



where i,j — L,R and are the chirality projection operators. Note that the Lagrangian of 
eqn.(l) is by no means a comprehensive one and similar operators involving the quarks alone 
(or the leptons alone) would also exist. However, for our purpose, it would suffice to consider 
only eqn.(l). Within this limited sphere of applicability, the strength of the interaction may be 
entirely absorbed in the scale A, and the couplings rjij and ^ij canonically normalized to ±1. 

Though the lagrangian presented above eqn.(l) is so-called 4- fermion contact interaction 
lagrangian but there are other theories which can give rise to such an effective interaction 
lagrangian. As is well known, a four-fermion process mediated by a particle with a mass signif- 
icantly higher than the energy transfer can be well approximated by a contact interaction [10] 
term with a generic form as in eqn.(l). For examples, theories with extended gauge sectors, 
leptoquarks [13], sfermion exchange in a supersymmetric theory with broken i?-parity [14] etc. 
are the theories which can give rise to such kind of effective interaction lagrangian by integrat- 
ing out fields with masses ^ A [15]. In those eqn.(l) arc just the lowest order (in A~^) ones 
among the series of such higher-dimensional terms. 

Such operators, in principle, could lead to significant phenomenological consequences in 
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collider experiments, whether e+e" [17,18], eP [16] or hadronic. Given the higher-dimensional 
nature of C, the fractional deviation over the SM expectations would be concentrated more at 
higher invariant masses. Some of the best constraints on compositeness, for example, came from 
the CDF [19] experiments and more recent measurement of the Drell-Yan cross section [20] at 
high invariant masses set the most stringent limits [21,22] on contact interactions of the type 
given in eqn.(l). For the e'^e~ collider, more recent constraints on compositeness, for example 
came from the OPAL [18] A ^ 1.6-3.4 TeV within the VA-iype interaction scenario. 

Recently, we have done the NLO QCD correction [23] in the context of hadron collider by 
taking into account of this effective lagrangian eqn.(l). In this article, we are going to do similar 
type of calculation in Linear Collider. In e^e~ annihilation, the perturbative QCD predicts 
only parton cross section but experimentally one measured only hadrons though hadronisation 
process known only phenomenologically. Because of this limited knowledge of hadronisation 
process one can not directly relate theory and experiment. Since we measured only hadrons in 
the final state one should include the higher order QCD corrections (which include the more 
partons in the final state) to the lowest order one to get the better result. In SM, people 
have done their calculation of next-to-leading order (NLO) [24,25] and next-to-next-to-leading 
order (NNLO) [26] QCD corrections for the dijet production in e^e~ annihilation. However, 
no calculations exist for the higher order QCD corrections to cross sections mediated by a 
generic contact interaction. Consequently, all extant collider studies of contact interaction have 
either been based on just the tree level calculations, or, in some cases, assumed that the higher 
order corrections are exactly the same as in the SM. Clearly, this is an unsatisfactory state of 
affairs and, in this paper, we aim to rectify this by calculating the next-to-leading order QCD 
corrections for both the FA-type and the 5'P-type contact interactions. 

Being nonrenormalisable of such theories, the current-current form of the lagrangian allow 
us to calculate reliable QCD corrections because of the fact that only one current consists of 
coloured field.This holds not only for the specific interaction in question, but also for other 
theories that satisfy the abovementioned criterion [28] . 

The rest of the article is organised as follows. In Section 2, we present our results for the 
LO and NLO cross section for resolved two parton case only (stated otherwise). The resolved 
three parton cross section will be divergence free and hence can be evaluated numerically in 
4-dimension. Here we also present only the total three parton cross section. In Section 3 we 
discuss the numerical results and finally, we summarise in Section 4. 



2 NLO corrections 

Before going to actual calculation of jet cross section, it is necessary to define cross section in 
perturbation theory. In perturbative QCD, each outgoing hard parton regarded as one jet that 
means one has to apply a jet resolution criterion to each outgoing partons to define a jet. In 
other words, a proper definition is to introduce a parton resolution criterion to define when a 
parton is resolved either as a single parton or as a cluster of partons. Consider the process e+e^ 
annihilate to quarks and gluon i.e. e+e~ — > q{pq) (l{Pq) g{Pg)- This process can be thought of 
as lowest order three jet production or higher order dijet production depending upon how we 
define the jet resolution criterion. One possible jet definition is a minimum mass cut so that the 
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invariant mass of pair of jet must be larger than experimentally defined one Smin- In the above 
mentioned process, there are only three invariant mass = (pi +Pj)^,i,j — Q,q,g- Therefore 
the lowest order three jet cross section is defined by 

da{e'^e'' — > 3 jets) — Qda{e'^e'' — > qqg) (2) 

and is the jet resolution criterion for the three jet final state defined by 

where 9{x) = 1 for a; > and otherwise. For the dijet cross section, it will represent the 
next-to-leading order (NLO) dijet cross section (where one of the partons (gluon, say) is either 
soft or collinear to other partons) and the jet resolution criterion for the two-jet final state 
would be the only one of the invariant masses {sqq say) larger than the experimentally defining 
cut Srfjiiji i.e. 

da{e^e^ 2 jets) = Qda{e^e^ qqg) (4) 

where 

Q ^ 0{Smin- Sqq). (5) 

We consider the process e'^e~ annihilation into a quark- antiquark pair i.e. e'^e~ qq in 
the context of generic contact interaction as defined by the lagrangian eqn.(l). In the presence 
of scalar-pseudoscalar (SP) type contact interaction, the leading order differential cross section 
for the above process is given by 

dcosd ~ 2 2^i,j=L,R\'itj\ \4, V / 

and for the vector-axial-vector (VA) type contact interaction, the leading order differential 
cross section for that same process (stated above) will be same as the standard model one and 
is given by (for completeness) 



d cos 2 s 



{\fLL\' + \fRR\') U] +{\fLR\'+\fRL\ 
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where t = — 1(1 — cos^), u — — 1(1 + cos^) and a is the electromagnetic couphng constant. 
fijihj ^ L, R) are given by 

fij = QiQq + gigjxis) + v^j^■ 
{8) 

X{s) = s/{s-Ml + tMzTz) 
The left-handed and right handed couplings g{^ and g{i of the fermion to Z-boson are given by, 

9l = . n - Qf Sin' ^w) , g{i = . n ( ' Qf ^w) (9) 

sm 9w cos 9w^ ' sm d^^ cos 6'^^^ ^ ' 

where e is the electron charge, Q/ is the electric charge in units of |e| of the fermion /, // is the 
third component of weakisospin and Q-^ is the electroweak mixing angle. 

At the leading order (LO) dijet calculation is much more simpler (as calculated above 
eqns.(6,7) than the next-to-leading order. At NLO, it requires careful treatment of cancellation 
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of divergences (soft and coUinear divergences against the divergences stemming from virtual 
corrections). The divergences coming from the fact that at NLO a parton can only be defined 
through a resolution criterion. There are many form of this resolution criterion. We have used 
the invariant mass resolution criterion. That is, if the invariant mass of the two parton less 
than the invariant mass resolution {smin) then these two parton considered as unresolved parton 
and treated as a one parton (or one jet) by integrating out the unresolved phase space which 
separates the soft and coUinear region of phase space from the resolved bremsstrahlung phase 
space. After adding this unresolved soft and coUinear contribution with virtual corrections, it 
becomes finite. 

For dijet production, the order as correction receive two contributions, one is from resolved 
two parton cross section (which is purely virtual contribution) and other one is the lowest order 
unresolved three parton cross section (calculated in the soft and coUinear limit) where one 
soft and/or coUinear parton clusters with one hard parton to form one jet. These soft and/or 
coUinear divergences can be isolated and it is easy to show that these divergences analytically 
cancel with soft and coUinear singularities coming from virtual correction of resolved two part 
process upto order Smin (or ymin = Smin/sij) where virtual gluon in the loop becomes soft. The 
lowest order unresolved three parton cross section will be leading order cross section multiplied 
by some functions F^'^^ (a function of all the singularities (soft and coUinear), Smin (or Umin) 
and the factorisation scale (a*f))- This F^'^'^ is same for type or 5'P-type theory and can 
be found out in the literature (for example see [24]). The lowest order three parton cross section 
is given by. 



"''^h ^ 27rr(l-e) 
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(10) 



(77 = SP, VA) 



The virtual corrections to the resolved two parton process is also available in the literature 
(for example they can be read from Refs. [27] for S'P-type and Refs. [24,28] for VA type theory). 
We follow the dimensional regularisation procedure to regulate all the divergences in d = 4 — 2e 
dimensions and we use the MS scheme to remove the ultra-violet divergence. For completeness, 
they are given by 



where 



VA 



r(l-e) \A%IJ? 



(77 = SP, VA) 



-8 + 7r' 



(11) 
(12) 



SP 



r(l-e) \47rn'^ 



2 + 7r' 



(13) 



From eqns.(10,ll), it is clear that the left over soft and/or coUinear divergences (from the 
virtual corrections) cancel against the soft and coUinear divergences coming from unresolved 
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three parton process. Therefore next-to-leading order cross section for exclusive dijet production 
are given by 



da 



(0) 
VA 



X 



1 + 



27rr(l-e) 



^ - 1 - 3 In (yrnin) " 2 In^ (ymin) I + (^(z/min) 



(14) 



and 



daPp X 



1 + 



27rr(l-e) I 3 



+ 5 - 3 In (ymin) - 2 In^ (ymin) ) +0(y„ 



(15) 



which are finite as e — > 0. Here we have considered fx — ixp- By integrating over resolved three 
particle phase space, one can get the 0{ymin) correction of the above eqns.(15,14). The resolved 
three parton cross section for the VA-tjpe contact interaction will be same as SM which is 
available in literature (see for example Ref. [26,30]) and is given by. 



a 



VA 
3- jet 



(0) a^Cp 
^VA 27rr(l-e) 



- ^ + I + 3 In [yrain) + 2 In^ [y^^) + 



(16) 



For the S'P-type contact interaction we have calculated the resolved three parton cross section 
as given below 



SP 
'3- jet 



— a 



(0) a.Cp 



SP 27rr(l-e) 



^ + I + 3 In (ymin) + 2 In^ (jjmin) + 0(ymin) 



(17) 



From eqns.(14,16) (or eqns.(15,17)), it is clear that the theoretical results for resolved two parton 
and three parton depend strongly on an arbitrary parameter ymin- Any physical observable 
should not depend on this arbitrary parameter. However for physical 2-jet NLO cross section, 
both two parton and three parton cross section will contribute and hence it is independent 
of this arbitrariness. This also ensures the KLN (Kinoshita-Lee-Nauenberg) theorem that the 
fully inclusive e+e" cross section is finite as quark mass goes to zero (i.e. free of mass singular). 



3 Results and Discussions 



In this section, we present numerical result for ILC. We choose the contact interaction scale (A) 
to be 2 TeV and the center of mass energy to be = 500 GeV. As is well known, the higher 
order QCD correction reduces the uncertainties related to scale choice (the renormahsation scale 
(//) and the factorisation scale (/^f)). For the NLO jet calculation, the analytical result does 
not depend on any of these scale explicitly [29] . The scale dependence comes through the strong 
coupling constant ^^(/i^). We have used the NLO as(/i^) for the NLO analysis and the scale we 
choose, both renormahsation and factorisation scale to be ii (iip) = Pt (otherwise stated). We 
have also shown that both two-parton and three parton result strongly depends on the cut-off 
ymin- For very small values of ymin-i two parton cross section become negative and three parton 
cross section become large positive (because of these terms \n{ymin)^ lii^(ymm))- For large 
enough ymin: both the parts produce the meaningful result. Our analytical result is valid for 
small ymin region since we have neglected the term 0{ymin) in the integration. Therefore ymin 
should be much less than one {ymin « 1)- For this reason, we choose ymin — 0.01 (detailed 
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discussion can be found in the literature [24,26]) for all the differential distributions and the 
total cross section. Furthermore, in presenting our results, we shall consider only one of the 
couphngs TjAB and ^ab to be non-zero and of unit strength. 

For the sake of convenience, we parametrize the cross section as 

c" = C"SM + c"intf + c"r72 (for the VA case) 

(18) 

a = (JsM + ^"^2 (for the SP case) 

and similarly for the differential cross sections. This has the advantage in that the total cross 
sections, for an arbitrary value of A can be easily reconstructed. We also take care of the 
so-called initial state radiation(lSR) effect [31] through out our numerical analysis otherwise 
stated. 
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Figure 1: The dijet production cross section as a function of invariant mass Mjj for y/S = 500 
GeV and A = 2 TeV. The scale is chosen to be /i(/iF) = Mjj. The upper set of lines (blue and 
pink) are for SP-type pure contact interaction (right most, lower panel). 

In figure 1 we have plotted the total cross section as a function of invariant mass of dijet 
(Mjj). For dijet production, the dijet invariant mass is same as the effective center of mass 
energy s (due to the ISR effect). The rapid fall of asM with Mjj is due to where as the 
interference cross section {(Jintf) is independent of s (eqn.(7)) and hence is almost constant. The 
slow fall of (Jintf reflects the higher dimensional nature of the contact interaction lagrangian. 
Here we have chosen the scale to be ^ {^p) = Mjj. The pure contact interaction cross section 
(7^2(^2) increases with Mjj due to fact that it is proportional to the s (eqns.(7,6)). The yA-type 
contact interaction cross section increases very slowly with Mjj because of its V — A current 
structure compare to S'P-type contact interaction cross section. Consequently larger value of 
Mjj, contact interaction dominates over the SM one. From figure 1, it is clear that the cross 
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section is flavor dependent whereas the purely contact cross section o"^2(^2) is flavor independent 
as we expected. The same argument holds for rest of the analysis. 

Figure 2 shows the angular distribution between beam axis and the jet axis. This distri- 
bution is almost constant for S'P-type contact interaction because of the fact that the leading 
order cross section is independent of 6 (which is typical characteristic of scalar vertex) as we 
see in eqn.(6). This argument holds not only for leading order result but for higher order 
corrections result as well. The small variation in the NLO result (yellow line) is due to scale 
variation through the ^^(/i^). Whereas this is not so for V^A-type contact interaction. V A-iy^ie 
contact interaction does depend on the 9. These differential distributions for both l^A-type 
and S'P-type contact interaction dominates over the SM piece. For the pure l^A-type contact 
interaction, the angular distributions are different for LL{RR) and LR{RL) sector. This is 
because of the sign of cos^ are different (eqn. (7)). In other words, their chirality structures are 
quite different. This is also true for Pr-distributions (see figure 3). 




cose cose cose 

Figure 2: The angular distribution for dijet production at ^/S = 500 GeV and A = 2 TeV. The 
scale fiifip) = Pt- Here fl is a scale factor. 

In figure 3 we have shown the transverse momentum (Pr)-distribution of a single jet. Since 
at the leading order, the transverse momentum of the two jets balanced each other (P^^ = 
= Pt)- For the NLO, the unobserved third parton can be taken infinitely soft for IR safety 
which is the artifact of fixed order perturbation theory. Therefore this momentum relation 
still holds even at NLO. From the figure it is clear that as Pt approaches towards \fS/2 
the differential cross section approaches infinity as we expected. Though the interference Pt 
differential distribution of RL and LR are of the same order or less but LL and RR distributions 
are of the order one more than the SM whereas for pure contact interaction it is ~ 10^ over the 
SM for both the cases {SP as well as VA). 
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Figure 3: Single jet Pt differential distribution for \fS = 500 GeF, A = 2 TeV and Umin = 0.01. 
The scale fiifip) = Pt- ^ is a scale factor. 



In figure 4 , we have plotted dijet cross section (only 0{as)) and resolved 3-jet cross section 
as function of i/min (without ISR effect). From the figure, one can easily see that large values of 
ymin, both two-parton and three parton cross section produce the meaningful result compared 
to very small valued i/min- We have also checked numerically that the sum of these dijet 
and resolved three jet cross section is independent of i/min (as it is cleared from the analytic 
structure in Section 2) which essentially reproduce the inclusive results for e~^e~ annihilation 
to quark-antiquark pair upto 0{as)- 



4 Conclusion 

To summarise, we have performed a systematic calculation of the next-to-leading order QCD 
corrections for the dijet production via e~^e~ annihilation in theories with contact interactions. 
Contrary to naive expectations, we demonstrate explicitly that the QCD corrections are mean- 
ingful and reliable in the sense that no undetermined parameters need be introduced. 

For the VA-tjpe interactions, the analytical structure of the corrections are similar to those 
for the SM. However, a significant dependence on the fiavor structure is found and needs to 
be carefully accounted for in obtaining any experimental bounds. For the SP-tjpe interaction, 
not only are the analytical results quite different, but the results are typically larger than those 
within the SM. Finally, we have investigated the sensitivity of our results to i/min (or s^m the 
invariant mass cut). 
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Figure 4: Dijet and 3-jet cross section as a function of Umin for y/s = 500 GeV, A = 2 TeV. 
The scale f^ifip) = Pt- ^ is a scale factor. Here "intf" signifies the interference cross section 
and "cont" signifies the pure contact interaction cross section. 
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